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Abstract 



^ . The idea of "dynamically" generated parton distribution functions, based on regular initial con- 
ditions at low momentum scale, is reanalyzed with particular emphasize paid to its compatibility 

' with the factorization mechanism. Basic consequences of this approach are discussed and compared 

CN . to those of the conventional approach, employing singular initial distribution functions. 

T— I ■ 

^ ■ 1 Introduction 
(N 

' One of the highlights of the recent International Conference on High Energy Physics in Brussels 

■ has been the remarkable success of "dynamically" generated parton distribution functions (DGPD) 

\^ i advocated by Gliick, Reya and Vogt (GRV) in predicting the rapid rise of proton structure 

0> ; function ^2^(2;, Q) at low x, observed at HERA §,0. Q, 

r"| ' The GRV group is one of three main groups (the other two being the Durham (MRS) lll|-|3| 



and CTEQ |T^,[T^ ones), which systematically analyze hard scattering data within the framework 
Q^l of perturbative QCD. What distinguishes GRV approach from those of the other two groups is their 

claim that the DGPD are more than just parameterizations of our inability to compute structure 
functions directly from first principles. GRV argue that by imposing certain condition on the 
initial parton distributions at low momentum scale, one obtains more predictive results. Without 
^ ■ this additional theoretical input the conventional parameterizations, using a moderate initial scale 

^ . Qo ^ 2 GeV, are unstable when extrapolated to low x region. For that reason both the MRS and 

CTEQ groups usually present several sets of such parameterizations, differing just in low x region. 

The idea of DGPD is intuitively appealing and actually almost as old as QCD itself [^. Con- 
fronted with growing amount and variety of data, it has, however, undergone significant modifica- 
tions 0,^,0] and in the process lost most of its original appeal. As the GRV approach relies on very 
low initial scale in the range 0.5 — 0.6 GeV, it has been met with reservations and scepticism p^l8|. 



In response to this criticism and in order to bring further arguments in favour of their approach, 
GRV have included in their recent paper |^ an extensive discussion of several of these points. 

To relate physics of short distances, the true realm of perturbative QCD, to that of distances 
comparable to the proton size would certainly represent a major achievement. The purpose of this 
paper is to discuss whether this can really be done in the way suggested in ||3|-|8|. Throughout this 
paper I shall concentrate on the analysis of the basic assumptions and consequences of the GRV 
approach, with only occasional reference to comparison with experimental data. 

The paper is organized as follows. In the next section I shall briefly recall the development 
of the idea of DGPD, from its inception |§ up to the present status 0]. In Section 3 the ap- 
plicability of perturbative QCD at distances as large as 0.4 fm will be discussed. In particular 

^Throughout the paper the term "distribution" stands for distribution function. 
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1 shall comment on the implications and interpretation of recent lattice calculations [|T9[, quoted 
in 0. The indispensable role of power corrections in going from short distances (where partons 
live) to distances comparable to the proton size (where the appropriate degrees of freedom are the 
constituent quarks) is emphasized in Section 4. In Section 5 the compatibility of the DGPD with 
the factorization mechanism is discussed in detail. In particular it is shown why it is very difficult 
for gluons and sea quarks to be valence-like. This discussion also shows how the conventional 
parameterizations, based on singular input distributions, avoid this problem. In Section 6 results 
of the conventional approach in the small x region are briefly reviewed and cast into a simple form 
suitable for the comparison with the DGPD. This comparison, carried out in Section 7, identifies 
two basic signatures of the DGPD. Throughout the paper I adopt the notation in which the QCD 
couplant a = ^ satisfies the usual RG equation 

da(M, RS) ^ _^^2(^^ RS) (l + ca(M, RS) + C2a\M, RS) + • • •) , (1) 
d m M V / 

where b, c are the first two, universal, coefficients ^ 

, 33 -2n/- 153- 19n. 

b = ^, c= ^, 2) 

6 ' 66-4n/ ' ^ ^ 

while all the higher order coefficients ; A; > 2 in (|l|) are free parameters, defining the so called 
renormalization convention (RC) |2^. Together with the specification of the initial condition on 
the solution of (|l|) they define the renormalization scheme (RS). 

2 The evolution of the idea of dynamically generated partons 

The original idea of ^ was to generate parton distributions at large momentum scales, where 
experimental data are available, by means of the DGLAP leading order |^, leading twist evolution 
equations, starting at some small momentum scale /i ~ 0.55 GeV from purely valence-like quark 
distributions, with vanishing light sea and gluon ones Q 

G{x, /i) = u{x, n) = d{x, fj,) = s(x, n) = s(x, ^) = 0. (3) 

The quark distributions at the initial scale /i, obtained by backward evolution from measured 
structure function (x, Q) at Qo — 3 GeV^, were constrained to satisfy quark number sum rule 

1 

dx [n(x, /i) + ci(x, /i)] = 3, (4) 





which provides a fundamental bridge between the parton model of Feynman and the old nonrel- 
ativistic "quasinuclear colored model" of Gell-Mann, Zweig, Greenberg, Lipkin and others. The 
scale /i was fixed by imposing the momentum sum rule 

1 

dxx [n(x, /x) + (i(x, /x)] = 1. (5) 





^Assuming QCD with 3 colors and n/ massless quark flavors. 



^In later GRV papers also the NLO DGLAP evolution equations were used. 
*In the rest of this paper Qo is used for the general initial scale, while the symbol ^ is reserved for the initial scale 
within the GRV approach, i.e. the one at which the parton distributions become valence-like. 

^In GRV approach heavy quarks c, h and t are not considered as intrinsic partons in the nucleon, but are produced 
from intrinsic gluons via the boson-gluon fusion mechanism Q . 
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In 1976 there were too few data on hard scattering processes to test the DGPD thoroughly. 
With more and increasingly accurate data becoming available in late eighties, the GRV were forced 
to modify their original idea by allowing nonvanishing valence-like gluon distribution fi) as 
well Moreover, G{x,fi) was related to the input valence-like quark distributions as follows: 

HQ 

G{x,n) = — [u{x,n) +d{x,n)], q{x,n) = 0.. (6) 

In /i and no were fixed by means of the momentum sum rule at the initial scale /i, now including 
also the gluon contribution, together with the comparison of theoretical predictions with data on 
direct photon production. They got about the same as in Q and nc = 2. 

Confronted with still more data GRV had finally to include in their initial parton distributions 
also the valence-like nonstrange sea Q so that the momentum sum rule now reads 

1 _ -, 

dxx Uy{x, fi) + dv{x, fj,) + 2u{x, fi) + 2d{x, fj,) + G{x, fj,) =1. (7) 

In one of their latest NLO global analysis, published in ||7|, /i = 0.58 GeV and 

xu^{x,^i) = 0.988x°-^''^ (l + 1.58V^+ 2.58X + 18.1x^/2) (l-x)^■^^ (8) 
xdy{x,^i) = 0.182x°-^^^ (l + 2.5lV^+ 25.0X + 11.4x^/2) (l-x)^■^^^ (9) 

x{u + d){x,fi) = 1.09x°-3(l + 2.65x)(l -x)^■^^ (10) 
xG(x,^) = 26.2x^-9(1 -x)^-°, (11) 
xs{x,fi) = xs(x,/x) = 0. (12) 

3 Does perturbative QCD make sense at 0.4 fermi? 

The initial scale fi = 0.58 GeV corresponds to a distance 0.37 fermi. As pointed out by a number 
of authors |1^18| such distances are probably too large for a meaningful purely perturbative 



treatment. It is fair to say that GRV do not trust their results at such low momentum scales 
but claim they are good approximations only above somewhat higher scale fipevt = 0.75 GeV ||7|. 
However, even the latter value seems too low for the applicability of leading twist, low order (LO 
or NLO) perturbative QCD. The following paragraphs are intended to throw some light on this 
problem. 

In this context let me first comment on the results of recent lattice calculation |l^ , quoted in ||7| . 
According to Q these results "confirm the perturbative NLO (2 loop) predictions for as{Q) down to 
Q = 0.55 GeV." This claim relies on Fig. 1, taken from jl^, where the results of nonperturbative 
lattice evaluation f\ of the QCD running coupling a!f**(g) is plotted as a function of q in the 
region q G (0.5,14) GeV. The agreement between a^^**(g) and the curve corresponding to 2-loop 
perturbative /3-function, i.e. the solution of (||) with only the first two universal terms on its r.h.s., 
down to 0.5 GeV is indeed remarkable. However, what is demonstrated by Fig. 1 is merely the 
fact that a particularly defined lattice couplant a!f** coincides with the couplant defined in the so 
called 't Hooft RC |^ down to q = 0.55 GeV. It is worth emphasizing that even on the lattice there 
is no unique "nonperturbative" /3-function and, consequently, no unique nonperturbative couplant 
Q,iatt| \Y]2iig asymptotic freedom of QCD guarantees that at short distances couplants in different 
RS coincide, they may be arbitrarily far apart at large ones. Fig. 1 contains an interesting evidence 

®The definition of a'/" used in is rattier involved and is tfierefore not mentioned liere. 
^In tliis RC all nonunique /3-function coefficients Cj,j > 2 are set to zero by definition. 
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for the closeness of two (out of an infinite number) definitions of the couplant, but tells us nothing 
about the applicability of perturbation theory in any of them. The authors of |jl^ are well aware of 
this limitation and on page 495 of [jl^ therefore write: "We would like to emphasize, however, that 
our results do not prove that perturbation theory provides a good approximation to all quantities 
of interest up to couplings as large as 3.48. Such a general statement is bound to be false and the 
running coupling in our scheme may very well turn to be an exceptional case." 

To illustrate the importance of nonperturbative effects at distances 0.3—0.4 fermi, let us consider 
the magnitude F{r) of the force between two static quarks at a distance r. This quantity has been 
extensively studied on the lattice and is usually written as the sum 



^(^) = o + ^ = ^P(^) + ^np, (13) 



4 (r) 

where the first term, dominant at short distances, comes purely from perturbation theory while the 
second describes the nonperturbative, long range confining force with k denoting the string tension. 
The couplant af^ in the numerator of (ll^) is related to a^^^^, mentioned above, as follows 



af(M) = a'/"(M)(l + A;ia!f"(M) + ---) , fci = 1.33776. (14) 



Evaluating both terms in ([ij) for k = (0.48 GeV)^ and two values of M (M = /x = 0.55 GeV, 
corresponding to r = 0.37 fm and M = /ipcrt = 0.75 GeV, corresponding to rpert = 0.27 fm), using 
the values of a^^^^ from [ff9| , we find 



Fp(0.37 fm) = 0.17 GeV^, Fp{0.27 fm) = 0.27 GeV^, Fnp = 0.24 GeV^. (15) 

This is the kind of comparison which really tells us how important is the perturbative contribution 
to a particular physical quantity, in this case the interquark force ([T^). For this quantity the 
nonperturbative contribution clearly dominates over the perturbative one at the distance r = 
0.37 fm and is roughly equal to it at r = 0.27 fm, the distance at which dynamical perturbative 
predictions should, according to [^, "become reliable and experimentally relevant". The above 
example is merely an illustration and the relative importance of perturbative and nonperturbative 
parts may well depend on the physical quantity in question, but it at least gives some indication that 
at //pert = 0.75 GeV perturbation contributions can hardly be expected to be a good approximation 
to the full results. 

One of the basic features of perturbation theory at low momentum scales (large distances) is the 
increasing sensitivity of finite order perturbative approximants to the choice of the renormalization 
and factorization schemes and renormalization and factorization scales |^, which makes perturbative 
predictions progressively more ambiguous in this region. To illustrate the crucial importance at 
large distances of higher order terms in purely perturbative expansions let me briefly recall the 
essence of Ref. There the familiar i2-ratio in e"'"e~ annihilation into hadrons at the center of 
mass energy Q 

He..- (9) ^ '^fiy ' = 3 f f e^) (1 + r(Q)). (16) 



(T(e+e — > /i+/x 
where in perturbative QCD 

r{Q) = a(M, RS) [l + ri(Q/M, RS)a(M, RS) + r2(Q/M, RS)a2(M, RS) + • • •] (17) 

is investigated in the infrared region. Their analysis has two important ingredients 
*The value 3.48 quoted above corresponds to g^. 

^In this note these two in principle diflterent scales will be identified and the dependence on the choice of factor- 
ization scheme disregarded. 
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• the use of the NNLO approximation to (|T7|), with C2 chosen by means of the PMS |2C], 

• smearing of Re+e-{Q) over some interval A of Q: 

«,..-«3.A)4/;d.^^^^^ (18) 

The second ingredient of this procedure is vital as the detailed structure of Re+e- (Q) in the res- 
onance regions is clearly beyond the reach of perturbative QCD. Nevertheless the fact that after 
the smearing and for not too small values of A, the agreement between (|18|) and the data is quite 



good, as documented by Fig. 6 of 1 21], is remarkable. This agreement depends crucially on the fact 
that C2, as chosen by the PMS, is negative, since for positive C2 the NNLO perturbative expansion 
for (|T^ ) blows up in the IR region, as does the NLO one. Also the magnitude of C2 is important 
as it determines the magnitude of R in the IR region. The success of such a procedure might look 
suspicious as in QCD sum rule approach resonance parameters (and thus also their contribution 
to the smeared spectra), are dual to power corrections. The observed agreement between data and 
NNLO perturbative approximation in the PMS approach can be interpreted as a signal that by 
an appropriate choice of C2 (or in general of the RC), it may be possible to include in some sense 
also the effects of these power corrections. Such an interplay of perturbative and nonperturbative 
contributions is quite plausible, as they actually coexist within the OPE. I have mentioned the 
analysis of |^l[ merely to emphasize that at large distances the inclusion of NNLO perturbative 
terms is probably indispensable for meaningful and reasonably complete description of physical 
quantities. 

The message of the previous paragraph may also have some relevance for "perturbative stability" 
observed within the GRV approach in the comparisons of the LO and NLO approximations to 
the leading twist DGLAP equations. Since the LO approximation cannot be associated with any 
well-defined RS, and the importance of higher order terms depends sensitively on the RS, such 
a comparison makes little quantitative sense. Only by comparing the NLO and the NNLO can 
such an information be obtained. Unfortunately, there are only a few simple quantities for which 
the NNLO calculations are available. One of them is the perturbative part (17) of the ratio (^). 
Evaluating just for illustration the first three known terms of (|17] ) at Upert — 

0.75 GeV, for three 

flavors and in the conventional MS RS, we find that they are roughly in the ratio 1 : 0.22 : 0.33! 
Not only is there no sign of perturbative stability for this quantity, but at such low scales the 
complicated problem of the presumable divergence of perturbation expansions in fixed RS becomes 
of utmost phenomenological importance. 



4 The interpretation of input parton distributions 

Let us now turn to the interpretation of parton distributions at the initial scale scale fi. While in Q 
they were considered to correspond to three constituent quarks, according to the latest paper Q 
the initial valence-like quark and gluon distributions "should rather be identified with current 
quark content of hadrons." This slight but crucial shift of interpretation should justify why the 
sea and gluon distributions do not vanish at the initial scale /i, as would seem appropriate for real 
constituents of the proton. Such an extension would still be reasonable were the additional sea 
quark and gluon valence-like initial distributions small admixtures to the basically three valence 
quark component of the nucleon. 

In Q the fact that two valence "constituent" gluons at /i ~ 0.5 GeV were required by the data 
was considered "fine" as "they may combine to give color and spin singlets as is required for the 
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nucleoli.". However, the success of the conventional SU(6) quark model relies on the fact that all 
color singlet combinations of three constituent quarks do exist in the nature, not only some of 
them! But in the color singlet state of three quarks and two gluons the latter do not have to couple 
to a color singlet. The system of three quarks and two gluons would have much richer spectrum 
of low-lying states than the state of mere three quarks. It would be a kind of "hybrid" states, 
suggested in the early eighties, but never found. In order to avoid these problems arguments would 
have to be invented to show why these two constituent gluons must (and not only may) couple to 
a color and spin singlet. I am not aware of any such argument. 

Similar problems arise for the initial distributions in summarized at the end of Section 2. 
Integrating over the initial distributions without the prefactor x to get the probabilities, we find 

J G{x, fj,)dx = 1, J (u{x, n) + d{x, iJ:)^ dx = 1.6. (19) 

The more accurate and copious data used in thus lead to the result that the initial parton 
distributions describe a system composed of 2.72 u quarks, 1.88 d quarks, 0.72 u antiquark, 0.88 
d antiquark and about one valence gluon. So again valence sea (anti)quarks and gluons are by no 
means a small admixture but on the contrary provide a dominant component of the initial parton 
distributions. 

According to GRV this is no cause for concern as perturbation theory is not expected to hold 
at the initial scale fj, « 0.55 GeV, but only above /Upert = 0.75 GeV. If, however, initial parton 
distributions do not describe at least approximately physics at the scale fi, what justifies then 
the adjective "dynamical"? In particular why to impose the fundamental sum rules (^) and (^), 
or d^, upon which the GRV approach is based? For instance, if the initial parton distributions 
q{x, fj,),q{x, fi) and G{x,fi) are irrelevant for physics, why should they satisfy the momentum sum 
rule (0), which expresses the fact that quarks, antiquarks and gluons carry together the whole 
momentum of the proton? And why should there be just two u and one d quarks at the scale /i? 
The answer GRV offer exploits the invariance of these sum rules under the LO DGLAP evolution 
equations and the fact their they do hold at short distances. This reasoning goes, however, against 
the very spirit of the DGPD, which I see in the possibility to use some known features of physics 
at long distances in order to predict physics at short ones. Imposing restrictions on the initial 
distributions in the situation when the latter have no physical meaning is basically a mathematical 
game with little physical content. Nevertheless, such a game can have nontrivial consequences, 
which, if confirmed by data, would signal some interesting physics behind the GRV approach and 
would justify it a posteriori. 

In my view the relation between constituent quarks and partons cannot be described by leading 
twist DGLAP evolution equations, even if these were taken to all orders, As pointed out in 
Section 2, the distance ~ 0.4 fm, which corresponds to /i ~ 0.55, is not much smaller that the 
approximate size of three constituent quark in the proton. The whole point of introducing the 
concept of constituent quark is that it represents the effective degree of freedom appropriate for 
describing the proton at distances comparable to its size. Contrary to the current quarks, which 
are associated with quark fields entering directly the QCD lagrangian, constituent quarks have no 
such firm basis and are merely an intuitively introduced concept, a kind of quasiparticle, which 
is reasonably well-defined only in the nonrelativistic quark model! The phenomenological success 
of this model suggests that in low momentum transfer processes proton behaves approximately as 

^"There is an misleading claim in [^] that their input distribution functions imply that "proton consists dominantly 
of valence quarks and valence-like gluons, with only 10% qq excitations (sea quarks)." The mentioned 10% concerns 
the momentum fraction carrie d by sea quarks and antiquarks, not the probabilities themselves. 

^^For related discussion see |22|. 



6 



composed of three constituent quarks, each with a mass of about 300 MeV. Constructing the fields 
corresponding to constituent quarks from those of the current quarks could probably be compared 
to the Bogolubov transformation between electrons and fermionic quasiparticles in the BCS theory 
of superconductivity. We expect the transition from short distances, where current partons are 
the right degrees of freedom, to large ones, where constituent quarks are the appropriate effective 
degrees of freedom, to be smooth, but involve complicated multiparton effects. In the framework 
of OPE such effects are described by multiparton distributions, which naturally appear as part 
of power corrections [23|. We cannot hope to get constituent quarks at low scale from quarks, 
antiquarks and gluons at large scales merely by means of the leading twist DGLAP evolution 
equations. 



5 Valence— like initial distributions and factorization 

In this Section the compatibility of the GRV approach with the mechanism of factorization of 
parallel singularities will be discussed. In order to make the discussion as clear as possible and to 
concentrate on the essence of the problem, a number of simplifications will be made. 

First, we shall be primarily interested in the low x domain, roughly x < 10~^. Secondly, all the 
considerations will be done within the LO DGLAP equations. The inclusion of the NLO corrections 
is not essential for any of the points discussed below. In the LO approximation the basic quantity 
of interest, proton structure function can be expressed in term of elementary quark 

distributions as follows: 



Q2) + Q^)) + ^ + d(x, Q2) + six, Q2) + s{x, Q' 



(20) 



Assuming SU(3) symmetry of the proton sea the r.h.s. of (20) can be written as a combination of 
the valence Uv, dy and the common sea D = Usea = c^sea = -Ssea distributions 



(x, Q2) = lxu^(x, g2) + ^xd^ix, Q2) + ^xD{x, Q^), 

y y o 

or, alternatively, as a sum of separate contributions from u and d quarks 



Fr{x,Q' 



-XUr, 



y 



1 



(x,Q2)+ a;Z)H(2;,Q2) + ^^^(^^g2)^ ^^(d)^^^g2) 



In terms of conventional moments of various functions (distribution, branching, etc.) 



(21) 



(22) 



/(n,Q)^ / x"/(x,Q)dx 

JO 

the LO DGLAP evolution equation for the nonsinglet quark distribution reads 

dgNs(f^,Q) "s(Q)d(o)/ n / r>\ 
dlnQ =^-P.VHgNs(n,Q), 

while for the quark singlet and gluon distributions we have a system of coupled equations 



dG(n,Q) ^ a.(Q) 
dlnQ TT 

d(g(n,Q) + g(n,Q)) ^ «s(Q) 
dlnQ vr 



p(S(n)G(n,Q) +P^;^(n)(g(n,Q) +g(n,Q)) 



,(0) 



2nfPf^{n)G{n,Q) + Pj^f{n){q{n,Q) +q{n,Q)) 



(23) 

(24) 

(25) 
(26) 



7 



As in the small x region the evolution of the gluon distribution is driven by the branching G — > 
G + G, we shall furthermore drop the second term on the r.h.s. of (^). Moments of the gluon 
distribution satisfy then the same kind of differential equation as quark nonsinglet distribution: 



dG{n,Q) 
dlnQ 



'^P§lin)Gin,Q) 



(27) 



and therefore also the corresponding solutions have the same form 



qmin,Q) = Aisisin) 
G{n,Q) = Acin) 



ca{Q) 



1 + ca{Q) 

ca{Q) 
1 + ca{Q) 



(28) 
(29) 



where A-^s{n), Ac{n) are unique finite constants, determining the asymptotic behaviour of the 
moments 5ns('T'5 Q), G(n, Q) as Q — > oo. According to the factorization mechanism \24\ these 
constants contain all the information on long range properties of the nucleon, uncalculable in 
perturbation theory. They represent one way of specifying the boundary conditions on the solution 
of evolution equations (|2^ ) and (|27|). Another and almost, but not entirely, equivalent way follows 
from taking the ratio of ( |28| ) and (|2^) at two different scales. In this case boundary conditions are 
specified at finite initial Qq and we have 



qNs{n,Q) 
G(n,Q) 



qmin,Qo 
G{n,Qo) 



a{Qo) l + ca(Q) 
a{Q) l + ca(Q, 



0) 



a{Qo) l + ca{Q) 



a{Q) l + ca{Qo) 



(30) 
(31) 



Let me first discuss the compatibility of the initial valence-like distributions with the factorization 
on the simpler case of the gluon density G{x,Q). Note that although they are of the same form, 
there is a profound difference between the solutions ( 25 ) and (p^) for n = (the moment giving the 



integral over the parton distributions). This is due to the fact that while Pqq^ (0) = 0, Pqq{0) = +oo! 
The former is a consequence of quark number conservation in the q ^ q + G branching, while the 
latter comes from the j spectrum of soft gluons. Eq. ([2^ ) implies that the integral over the 
nonsinglet quark density, (?ns(0) Q) = ^NS is independent of Q and provided it is finite at some Q 
it stays so everywhere. For the gluons there are two possibilities: 

a) Ag{0) > and then G{x, Q) cannot be valence-like at any Q for which a > since 



>{o), 



G(0,Q) =Ag(0) 



1 + ca{Q) 



ca{Q) 



Ag(0)(+oo) =oo. 



(32) 



b) Ag{Q) = 0, in which case (I 



is ill-defined, but ( |3lD applied to n = can still be used 
a(Qo) l + ca{Q) 



G(0,Q) = G(0,Qo) 



a{Q) 1 + ca((5o 



(33) 



If we now impose valence-like behavior on the gluon distribution at the initial Qq = /i, i.e. 
assume finite G(0,/x), (53) implies, due to monotonous behavior of the square bracket as a 
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function of Q, that G(0, Q) is a discontinuous function of the factorization scale Q at Q = fi: 



G(0,Q) = +00, yQ>i^ 

G{0,ii) = const. >0 (34) 
G(0,Q) =0 VQ</i. 

Similar discontinuity at the initial scale /x appears also for sea quarks and antiquarks. 

In the case b), realized in the GRV approach, it is difficult to understand why the divergence of 
^gg(O)' which is a purely perturbative phenomenon, should be accompanied by the vanishing of 
nonperturbative quantity Ag{0). For instance, in 4 — e dimensions PQQ{n,e) is finite and there is 

no obvious reason to expect Ag{0, e) = 0. Sending e — > we understand why Pqq{0, e) —>■ 00, but 
why should simultaneously Ag{0, e) — > 0? Although the lowest twist contribution may not provide 
a reliable description of the proton at the initial scale we are not free to impose arbitrary 
constraints on its properties. Despite these reservations, the case b) is certainly mathematically 
interesting option and I shall therefore in the rest of this Section analyze some if its consequences, 
in particular for the behavior of ^2*^(2;, Q) as a function of Q. 

In the small x region, to which we restrict our attention, the LO expression for the distribution 
of sea quarks within a single quark which at the scale n is described by the initial distribution 
6{1 — x), has the form ||2^ 

I^o(x,Q) = -fc2CV<^, (35) 
X 6 

where for 3 colors and Uf flavors 

+ - "-V'-'SCHIM c,^l (36) 

and l2iv) is the modified Bessel function. For general initial distribution q{x,^) we have 

"1 dy 



Jx y 



if ftc=e-"C [^y,(yW-P (37) 



where now w = \/48C ln(y/x). Note that for fixed x and C ^ 0, corresponding to Q ^ v ^ 
and Dq{x,Q) behaves as 

Do{x,Q) = -^C2e, (38) 
x 6 

vanishing for C = 0) at Q = fi. However, the physically relevant case of fixed Q > fj, and x ^ 
corresponds to slightly more complicated limit f — > co. Eq. (35) then implies 



z,„(.,g^)4fc.cV"<«W-^^ (3,) 

where the terms depending on v{x) induce additional dependence on x. This modifies slightly the 
behavior of Do{x,Q) at small x, but does not change its nonintegrability due to the dominant ^ 
factor. Consequently, Dq(x,Q), considered as a function of Q does not vanish uniformly in the 
whole interval x £ (0, 1) when Q ^ ^\ The same holds for D{x, Q). This nonuniform convergence 
means that for Q arbitrarily close to the input Qo = fi, there is always a region of x close to 
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X = 0, where approximately D{x,Q) oc 1/x. And it is this region which causes the divergence of 
the integral 

/ dxD{x,Q) = oo (40) 
Jo 

for any Q > fj,. In the conventional approach with singular, i.e. nonintegrable, initial distributions, 
this discontinuity is absent. Using eq. ( ^7|) and assuming for small x quark initial distribution in 
the form q{x, fi) = Ax^'^, A > 1, we find that for small x and Q ^jl^ the integral 

[' Ayq{y,Q)h!^ (41) 
Jx 

behaves differently than for A < 1. For A > 1 the integrand of (^) is a nonintegrable function of x 
in the interval (0, 1), finite lower integration bound is therefore crucial and we get for any Q > fi 

[\yq[y,Q)^-^x^-^ =^ D{x,Q)^- dq{y,fi)^Ax-\ (42) 
Jx A - 1 X Jx 

The singular initial distribution overrides the radiation pattern characterizing the emissions from 
individual quarks and the radiated sea quark distribution D{x, Q) is therefore of the same form as 
the initial q{x, ji). This well-known feature of singular initial distributions implies that for low x the 
form of the x-dependence is essentially independent of Q and the initial scale plays no exceptional 
role. 



6 Conventional partons in small x region 

In this section results of the conventional analysis based on singular initial distributions will be cast 
into a simple form suitable for the comparison with GRV results. In order to check the very essence 
of the GRV approach only the results based on the original initial conditions (^ with vanishing 
antiquark and gluon initial distributions will be discussed. Moreover, I shall concentrate on the 
small X region, where the approximations of the previous section are expected to hold |25|. In the 
approximation of neglecting the second term in ( p5| ) H{x,M) = xG{x,M) satisfies the equation 



dH{x,M) 
dlnM 



zP!,%{z)) H{x/z,M), 



where 



1 



+ 



1 



+ z(l - z) + 



/33 - 2n 



/ 



X 



V 36 



lU(i-z) 



Assuming H{x, M) in a singular factorizable form 

H{x,M) =x-^F{M,X), A>0 
and substituting ( |45| ) into ( |43| ) we get 
dF(M,A) 



dlnM 



a{M)F{M, A) C(A) 



dzz'Pt,%iz) 



where 

7(°)(A) 



dzz'pi,%{z 



6 



A(A + l) + (A + 2)(A + 3) 3 



A(A + l)(A + 2)(A + 3) 



1 



fc=0 



{2 + k){X + 2 + k)' 



(43) 



(44) 



(45) 



(46) 



(47) 
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extends the definition of the LO anomalous dimension ■j^^\n) to noninteger positive values of A. 
Note that 7(''^(A) is a decreasing function of its argument and negative for A above Aq ~ 0.85. The 
second term in the brackets of ( |46| ) is proportional to and can therefore be neglected in the 
small X region. F{M, A) then satisfies the equation 



which has a simple solution 



^i^|^ = a(Mh(°)(A)F(M,A), (48) 



F(M,A) = yl(a(M))-^' (49) 



where A is an arbitrary overall normalization constant. 

In the small x region the comparison of the GRV and conventional results can be made partic- 
ularly transparent by translating ( ^9| ) into the corresponding expression for the proton structure 
function F2^{x, Q"^) by means of a LO Prytz's relation ||2^ 



^l^g — = Ka{Q)H{x,Q), K=^- (50) 

This formula was shown to be a good approximation for F2 in the low x region and should be 
sufficient for our purposes. Alternatively, we could solve the DGLAP evolution equations for the 
coupled quark singlet and gluon distributions exactly and then approximate these solutions by the 
power-like behavior, but the procedure based on the combination of ( ]49| ) and (^) is much simpler. 
Anticipating also in the factorizable form 

Ff(x,Q)=x-^Fr(Q,A), (51) 

and using ( |50|) in combination with the explicit result ( |4^ ) for F{M, A) we get 

dF|P(M,A ) _ ^,^^1^^ 
dlnM 

wherefrom 



The positivity of F2^{M,\) requires positive C(A), which in turn implies A < Aq. The measured 
behavior of F2^ in the small x region is well within this limit. 



7 DGPD vs. conventional partons — numerical comparison 

Can the DGPD be distinguished from the conventional parton parameterizations at all? In order to 
identify reasonably unambiguous signatures of DGPD I shall concentrate on its "orthodox" version. 

In the first kind of comparisons GRV results for F2^{x, Q'^) were obtained via (^) from u and d 
quark valence-like initial distributions (^)-(|^ 0. The valence distributions were evolved by means 
of the standard DGLAP evolution equations using the method of Jacobi polynomials [ 
the sea parts (|37|) was used. Three light quarks were taken into account in generating the sea. 
In order to facilitate the absolute comparison between the two approaches without reference to 
experimental data, the constants A and A in (^) were related in such a way that the conventional 
results coincide with the GRV ones for x = 10~^ and = 5 GeV^. As a result A becomes a 

^^No essential features of the following comparisons depend on this particular choice of the initial distributions. 
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function of A. The choice of the normahzation point is, of course, to a large extent arbitrary, but 
none of the basic messages of this section depends on it. In Fig. 2a results of the GRV approach 
are plotted for 10 values of = o.35, 0.4, 0.6, 1, 2, 5, 10, 20, 50, 100 GeV^, i.e. starting very 
close to the initial fj? = 0.34 GeV'^. The thick solid curve corresponds to the initial F2^{x,fi). In 
Fig. 2b the corresponding sea distributions, which in the low x region are well approximated by 
the power-like behavior of the form x~^, are plotted. Figs. 2a,b nicely illustrate the way -F|''(x, Q) 
approaches the input function F2^{x^^) as Q ^ In particular the nonuniform convergence, 
discussed in Section 5, is clearly visible. As Q ^ fi'^ the region of increasing moves steadily 
to the left, but regardless of how close Q is to such a region eventually appears, leading to the 
discontinuity of the integral over quark distributions at the initial scale /i. 

Fitting for each the total GRV sea contribution in the interval x G (10~^, 10~^) to the form 
A{Q)x~'*'^'^^ leads to characteristic dependence of the parameters A{Q) and A(Q) on Q, displayed 
in Figs. 3a, b by thick solid curves. The corresponding power-like fits are shown as dotted curves 
in Figs. 2a,b. As Q ^ the interval where they provide good approximation to full GRV results 
shifts systematically to smaller values of x, reflecting the shift to the left of the relative importance 
of the sea component. 

In Fig. 4 GRV results are compared with those of the conventional approach, shown as dashed 
lines, for six values of A and the same values of as in Fig. 2. The lowest, thick solid curves 
of the GRV approach correspond again to the initial distributions ®-® at /i^ = 0.34 GeV^ and 
have thus no analog among the dashed curves. Note that for smaller values of A not all curves of 
the conventional approach fit into the frame of the plots. To summarize the message of the Fig. 
4, F^{Q^ A), defined in (51), is plotted in Fig. 3a as a function of for six values of A together 



with the corresponding function A{Q) of the GRV approach. 

In order to investigate the sensitivity of the parameters X{Q) and A{Q) in the GRV approach 
to the choice of the initial distributions, I have repeated the calculations presented in Fig. 2a for 
the set of u and d quark initial distributions of the form 

xq^{x,^i) = Ax'^il- xf (54) 

for two values of /? = 3,4 and five values of a = 0.1,0.3,0.5,0.7,0.9. In all cases the overall 
normalization factor A was fixed to give two u and one d quark at the initial scale ^i. Results 
corresponding to /? = 3 and the five mentioned values of a are displayed in Fig. 5a-e. The general 
pattern remains the same as in Fig. 2a, but there is a marked difference between the dependence 
of the exponent \{Q) (curves included in Fig. 3b) and the normalization factor A{Q) (see Fig. 5f) 
on Q and /3 (not shown). While the exponent X{Q) depends on the values of both a, /3 very weakly, 
A{Q) depends on a strongly and (3 moderately. The insensitivity of the exponent \{Q) to details 
of the initial distributions can be understood by closer inspection of the formula (p7|). The figures 
3-5 suggest two distinct features of the GRV results: 

• The characteristic dependence of the exponent \{Q). In the conventional approach A is 
arbitrary Q^-independent number, while in GRV one it is an almost unique function of Q, 
which starts rapidly from zero at the initial scale but then progressively slows down. The 
most sensitive region is clearly that close to the initial but even at large the characteristic 
dependence of GRV results persists. In the region between = 2 and = 100 GeV^, 
where leading twist can perhaps be trusted and data from HERA are available, A((5^) varies 
slowly in the range (0.25,0.38). Unfortunately the current experimental error on A at HERA 
is too large for a reliable discrimination of this dependence from a constant one, characterizing 
the conventional approach. To pin down the dependence expected in the GRV approach 
reliably would require lowering the experimental error on X{Q) to about 0.03 in a broad range 
of Q2. 
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• The dependence of the normahzation factors A{ff), F^iQ^X) respectively. Here the 
differences are much more pronounced and increase as we go to both smaU and large values of 
. They depend also much more on the choice of the initial distributions. In GRV approach 
we cannot go below the initial scale, while in the conventional approach there is no such strict 
limitation in the small x region. Fig. 3a indicates that the region close to the initial scale is 
particularly sensitive to the differences between the two approaches. 

The above comparisons show that the "orthodox" version of GRV partons could be easily 
distinguished from the conventional parameterizations, based on the singular initial distributions. 
However, once valence-like initial gluons and antiquarks are added the difference in the behavior 
of normalization factors A(ff) and F^iC^^X) largely disappears. This is documented by the fact 
that the latest GRV as well as the conventional parameterizations can accommodate new HERA 
data ||9|,[l^ at low x and in a broad range equally well. On the other hand the characteristic 

dependence of the exponent A(Q) remains essentially unchanged as it is still given by the basic 
radiation pattern (^). Similarly even in the refined versions of DGPD F'^ix^Q) approaches 
the initial structure function F^{x^\x) in same nonuniform manner as discussed in Section 5 and 
displayed in Figures 2 and 5. The available experimental data either reach sufficiently low values 
of x 10"^ but stop at > 2 GeV^ (HI and ZEUS at HERA), or reach low only touch 

the crucial region of a; < 10"'^ (E665 at Fermilab). Despite its limited x range recent E665 data 



indicate deviation from the latest GRV parameterizations for x below 10~^ and < 1 GeV^ [29 
New quantitative tests are expected from upgrades of HI and ZEUS detectors, which will extend 
the region of accessible down to a fraction of GeV^ and could soon throw some light on the 
questions discussed in this paper. 



8 Summary and conclusions 

In this paper I have discussed the basic idea and consequences of the DGPD from physical as well 
as mathematical points of view. I have argued that it runs into problems if we attempt to interpret 
the properties of initial distributions in physical terms. The only way to avoid these problems, and 
the one adopted by GRV, is to assume that the initial scale \x lies outside the range of validity of 
leading twist perturbative QCD. All the peculiar properties discussed in Sections 4,5 are then of no 
concern, but at the same time the approach looses its physical justification and becomes primarily 
an exercise in mathematics. Investigating the consequences of different initial parton distributions 
on the solutions of LO/NLO DGLAP evolution equations is an interesting mathematical problem 
of its own and I have therefore devoted the second part of this paper to it. In particular I have 
looked for signatures in the behavior of F^ix^Q'^) that would provide clears signals that GRV 
dynamics is at work. The only, but on the other hand rather unique signature of this kind is the 
characteristic dependence of the exponent in (51). To best way of pinning down this 

dependence at HERA is to extend the range below 1 GeV^, where the variation of is 
strongest, and study in detail the variation of in the whole available range. But even 
then a significant increase in the precision of measuring A is necessary before a definite conclusion 
on the Q^-dependence of A can be drawn. 
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Figure captions 

Fig. 1: The results of nonperturbative evaluation of particularly defined QCD running coupling 
on the lattice. Taken from Ref. [19]. 

Fig. 2: a) F2'{x,Q'^) as a function of x for the initial distribution (|8|)-(^) and 10 values of 
defined in the text. In the low x region the curves are ordered from below in order of 
increasing Q^. The dotted lines correspond to the power-like fits described in the text and 
the thick solid curve describes the initial F2^{x,fj,). b) The same as in a), but for the total 
sea component only. 

Fig. 3: Q^-dependence of the normalization factor A{Q) (a) and the exponent X{Q) (b). In 
a) the thick solid curve is given by the power-like fit of the GRV results, while the dashed 
lines, describing F|''(Q,A) of the conventional approach, correspond to six values of A = 
0.05, 0.1, 0.2, 0.3, 0.4 and 0.5. In b) the thick solid curve describes again the GRV result for 
the initial distributions @-(^), while the other curves correspond (from above in decreasing 
order of a) to five initial distributions (|5^ ) with /3 = 3 and a = 0.1, 0.3, 0.5, 0.7, 0.9. In b) 
similar curves for /3 = 4 would be essentially indistinguishable from those in the figure. 

Fig. 4: The comparison of the GRV (solid curves) and conventional (dashed curves) results for 
F2'^{x, Q^) separately for six values of A in the latter approach. 

Fig. 5: The same as in Fig. 2a (a-e) and Fig. 3a (f), but for /3 = 3 and five values of a in (p^). 
In f) the curves are ordered from above in order of increasing a. 
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